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Dirty, Skewed, and Backwards: The Smectic A-C Phase Transition in Aerogel 
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We study the smectic AC transition in anisotropic and uniaxial disordered environments, e.g., 
aerogel with an external field. We find very strange behavior of translational correlations: the 
low-temperature, lower-symmetry Smectic C phase is less translationally ordered than the high- 
temperature, higher- symmetry Smectic A phase, with short-ranged and algebraic translational corre- 
lations, respectively. Specifically, the A and C phase belong to the quasi-long-ranged translationally 
ordered " XY Bragg glass " and short-ranged translationally ordered " m = 1 Bragg glass " phase, 
respectively. The AC phase transition itself belongs to a new universality class, whose fixed points 
and exponents we find in a d = 5 — e expansion. 

PACS numbers: 61.30.Dk, 64.60.Fr, 64.70.Md, 82.70.-y 



I. INTRODUCTION 

Of all randomly pinned elastic media P, 0, S 3 liquid 
crystals in aerogel exhibit a phenomenon unique to 
themselves: anomalous elasticity. That is the scalings of 
their elastic energies are changed radically ( specifically, 
by non-trivial power laws ). 

However, there has been no previous work on phase 
transitions in pinned liquid crystal systems. In this paper 
we remedy this by treating the smectic A to smectic C 
( hereafter AC ) transition in an anisotropic, uniaxial 
disordered environment. Such an environment could be 
realized, e.g., by applying an electric or magnetic field to 
a liquid crystal in isotropic aerogel jS], or by stretching 
the aerogel uniaxially before absorbing the liquid crystal. 
We will hereafter refer to the special uniaxial direction 
as " along the applied field " or " the z-axis " . 

The AC phase transition separates the two novel glassy 
phases discovered in reference |^ . The high temperature 
phase (T > Tac) is the glassy analog of the smectic 
A phase of the pure problem, in that the layer normals 
lie, on average, along the applied field. This " random 
field XY smectic Bragg glass " (XYBG) phase is in the 
universality class of the random field XY model. The low 
temperature phase is the glassy analog of the smectic C 
phase, in that the layers normals make an angle 0{T) with 
the applied field. The experimentally measurable " tilt 
angle " 9{T) is the magnitude of the order parameter for 
the transition. This phase is in the universality class of 
the " m = 1 smectic Bragg glass "{m ~ 1) phase studied 
in 0. 

We call both of these phases " glassy" because the ran- 
dom environment (i.e., the aerogel) destroys long-ranged 
translational order in both. The extent of this destruc- 
tion, however, differs greatly between the two phases. 
Strikingly, it is the Zoti»-temperature, higher-symmetry, 
TO = 1 phase that has less translational order. In the 
XYBG phase, translational correlations are " quasi-long- 
ranged ", i.e., they decay as power laws with distance. 
In the TO = 1 phase, these correlations are short-ranged. 
This leads to radically different X-ray scattering signa- 
tures in the two phases which we will now describe. 



In the XYBG phase, the X-ray scattering intensity 
I{q) diverges near the smectic Bragg peaks, which occur 
at (f = nqaz for all n integer, where qa — —, with a the 
smectic layer spacing. This divergence is a power-law: 
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where a is a non-universal constant of order 1 and q± = 
\q-- qzz\. Note only the first 2 peaks (n = 1 and n = 2) 
actually diverge. In contrast, in the " glassy C " or " 
TO = 1 BG " phase, the peaks in the X-ray scattering 
intensity are broad, with I{q) finite for all q. 

As T ^ Tac from above (i.e., on the A side), the 
sharp peaks disappear in an unusual way. The peaks 
look broad for q 's sufficiently far from the Bragg peak 
position nqoz, while for g' s sufficiently close to the peak, 
they diverge per eqn.Q). " Sufficiently close " means that 
both \q± \ ^ 5(7^ (n,T), and {q^ — nqo\ ^ Sq'^(T), where 
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with - exp{A\T -Tac\ ), where n IS a univer- 
sal exponent calculated below and A is a non-universal 
constant. These predictions are illustrated in figures 1 
and 2. The divergence of z implies that, as T ^ 
the algebraic " spikes " on top of the broad short-ranged 
peaks get narrower and less intense, vanishing completely 
at Tac- Lowering temperature further leads only to the 
broad peaks of the Smectic C phase. This entire scenario 
of sharp peaks at high temperature and broad peaks at 
low temperature is very counterintuitive, and unlike al- 
most every other translationally ordered system • 

Why is the lower-symmetry C phase less translation- 
ally ordered than higher-symmetry A phase? In fact, it 
is precisely the new broken symmetry of the smectic C 
phase-that is, the tilt of the layer normal-that causes this. 
This is because, while the energetically preferred layer 
normal N in the smectic A phase is unique-it must point 
along z-there are infinitely many energetically preferred 
orientations of N in the C phase: N can lie anywhere 
on a cone making an angle 6{t) with z. As a result, the 
elasticity of the Smectic C phase is softer than that of 
the A phase because rotating the layers in such a way 
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FIG. 1: The -dependence of the X-ray scattering intensity 
for gx = in the smectic A phase. In the C phase, the sharp, 
power law peaks disappear, leaving only the broad scattering. 
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FIG. 2: The g_L -dependence of the X-ray scattering intensity 
for Qz = qo in the smectic A phase. Again, the sharp peak 
vanishes in the C phase. 



that TV rotates around this cone costs no energy. This 
exact symmetry of the elastic energy of the smectic C 
phase means that the direction perpendicular to the z- 
< N > plane (where < iV > is the mean of N; i.e., the 
direction of spontaneous tilt) becomes " soft ", that is, 
an easy direction for layer displacements to vary in. Pre- 
cisely such softness occurs (for different reasons) in the 
" m — 1 smectic " studied in 6] and, indeed, the elastic 
Hamiltonian for the C phase we find is identical to that 
studied for "to — 1 smectic" in :6|. Thus, we can simply 
transcribe the results of to this problem. 

In particular, positional fluctuations u of the layers 
about their optimal (tilted) positions obey < |u(9)p > = 



TC{q)-' + {AMq^s + ^h{q)ql + A.ql)C{q)-\ with 
C{q) = Bql+'-f{q}q1+K{q)q'^, where <> denotes a ther- 
mal average, and the overbar denotes an average over 
disorder. The anomalous quantities As^hi<t), Kio) and 
7(<z) obey 
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T-T^cr^^^"^) and 
qs is the component of q 



where qii — Eq^ — Fq^ with E cx 

perpendicular to the z-N plane, qx is the component of 
q within the z-N plane orthogonal to z and the critical 
exponents y± and r]t^K,c will be given later. The uni- 
versal exponents in (|3J) are given by = 2 — C^"^''^ ), 
^2 = 2 — ^ with the ?7's obeying the exact scaling re- 
lations 1 -\- rjs — ^ -\- 2fiK and rjh — 2 + rjs — rj^ — fiK- 
Numerical values for the ?7's have been obtained "^l by 
comparing e-expansions based on two different analytical 
continuations of this model to higher dimensions, giving 
f}K = 0.50 ± 0.03, r^-y ^ 0.26 ± 0.12, rj, = 0.132 ± 0.002 
and rjh = 1.372 ± 0.12. These predi ctions can be test ed 
by hght scattering, which measures < Ni{q)Nj{~q) > = 
qtqf< >■ The form given for < |M(g)p > also 

implies short-ranged translational order and, hence, 
broad Bragg peaks, as described earlier. 

We have studied the AC transition in an e = 5 — d 
expansion, where d is the dimension of space, and find 
that there is a stable fixed point, implying a second-order 
phase transition with universal critical behavior. In par- 
ticular, the tilt angle e{T) obeys OiT) = A{Tac ~ Tf, 
where, to leading order ine, (3~^~j^-\-0 (e^) . 

The specific heat exponent a — —-^ -\- O [e^) and the 
susceptibility exponent 7 = 0. 

The order parameter A^_l for this transition is the pro- 
jection of the smectic layer normals N perpendicular to 
the applied field. Above Tag, real space correlations of 
N± decay rapidly with distance, with correlation lengths 
^2 and parallel and perpendicular to the field re- 
spectively. Both diverge as power laws in (T — Tag)- 
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We find 



_L,2 OC 

-^ = l + f§ + 0(e2). 

We also find that the system exhibits anomalous elas- 
ticity right at Tag as well. Specifically, we find that, right 
at Tag-, the smectic layer bend modulus K vanishes as 
q^O according to the scaling laws 



qz , qz~> q\_ 



where the anisotropy exponent C = 2 — ^ and rj^ — 



Ck^^ + 0(e3) with Ck = '° ^ -.00353. 

that Ck < 0, which implies that K vanishes as q 
The anisotropy exponent also obeys C = — . 



Note 
^ 0. 
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The disordering effect of the random aerogel matrix 
can be quantified by disorder variances and de- 
scribing tilt and compressive stresses respectively. These 
variances also become anomalous, obeying (i = t, c) 



A. (g)=q-"7^(:^± 



— Vi 

9± . 



qz < 
qz > qS 



(5) 



Unlike the similar problem of a smectic A in isotropic 
aerogel and no field, the smectic layer compression mod- 
ulus B remains finite as — > at Tac- 

The exponents rjt and rjc are given by rjt ~ CAf^+0{e^) 

with Ca = ^^^"4y~^ « .01386 and = 2 - f -I- 0{e^). 

For T bigger than T^c the disorder variance Af c {q, T) 
and layer bend modulus K {q, T) are given by their T = 
Tac forms equations Q and Isj if either q±S,± ^ 1 or 
qz^z 3> 1. Otherwise (i.e., if both q±^± ^ 1 and q^^z <C 
1), At,c oc ^"l- oc (T^Tac)"'^'"-^ and K oc ^I" oc 

The critical exponents obey exact scaling relations: 
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(6) 
(7) 

(8) 



= C'^± aiid C = 2 — Note that a does not obey 
hyperscaling, due to the strongly relevant disorder. 

All of these exponents can be measured experimen- 
tally. The specific heat can, of course, be measured by 
the usual thermodynamic measurements. The sponta- 
neous tilt angle can not be deduced from the position 
of the smectic C Bragg peak, since that peak is broad. 

Fortunately, an alternative measure of 9q can be de- 
duced from the dielectric or diamagnetic susceptibility 
tensors Xij and e^. In the A phase, one of the principal 
axes of both tensors is along the applied field. In the 
C phase, this axis rotates away from the applied field 
due to the tipping of the layers. This rotation angle is 
proportional to 9q. 

The order parameter correlation lengths can be mea- 
sured by light scattering, which probes fluctuations in 
both the dielectric(eij) and diamagnetic(xy) susceptibil- 
ity tensors. The full form of the light scattering is ex- 
tremely rich and complicated; we will defer a complete 
description of it to a future publication. Here we will re- 
strict ourselves to pointing out that for qz = 0, the light 
scattering intensity scales like g^»?K->)t-4 ^J^, is 

independent of q±^ for {^^^) ^ (7_l <C ^J^, and scales 
like ^ for q^ « where ^ ^r^""^"* as 

T T^fj. Thus light scattering data should easily al- 
low determination of S.j_(T) and the combination of ex- 
ponents 2riK — rjt- Fitting the T-dependence of to 
(T - Tac)""^ then determines i^j^. 

We now briefly sketch the derivation of our results. 
Our starting point is an elastic energy for the layer dis- 
placement field u, which is the only soft variable in the 



problem, since the applied field locks the nematic direc- 
tor n 5]. For the smectic AC transition in a pure (i.e., 
disorder-free) system, with an applied field freezing the 
director out, Grinstein and Pelcovits showed that the 
appropriate elastic energy is: 



Bpure — I d> r 
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This model is very similar to that for a smectic A in 
the absence of an external field. However, because the ro- 
tational symmetry is broken due to the external field, a 
new term |V^up, which hardens the directions orthogo- 
nal to z, is generated. Since |V^u| is proportional to the 
tilt angle of the smectic layers, the coefhcient Dq{T) is 
positive in the A phase (favoring alignment of the layer 
normal with the applied field), and negative in the C 
phase (favoring tilt of the layers). Hence, by continuity, 
at T = Tac, T)q{T) vanishes. In what follows, we will 
assume that Dq{T) oc T — Tac near Tac- 

The other terms in © are simply those of the elastic 
theory of a smectic A in zero field, with one crucial ex- 
ception: in a smectic in zero field, rotation in variance 
requires that g = w = B, while for the AC in a non-zero 
field problem, at all temperatures, even at T = Tac, 
where D —> and softness is recovered, g and w are still 
free, because rotation invariance is still broken. 

To include the effects of the quenched disorder of the 
aerogel, we add to the pure Hamiltonian ^ random fields 
coupling to u and its gradients: giving us 



H = H, 



pure r 

/ (fr\h{r)-Vu + Vp{u - (t>{r)) 



(10) 



where h{r) is a quenched random field that for simplic- 
ity we take to be Gaussian zero distributed mean, and 
characterized by short-ranged anisotropic correlations: 



K (r-) h, (fO = [AtSi, + AcSt,] 6" (r - /) . (11) 

The field 0(r) is also a quenched random field with only 
short-ranged correlations, and is uniformly distributed 
between and a, the smectic layer spacing. The function 
Vp{u ~ (p) is periodic with period a. 

The physical interpretation of the quenched random 
fields h{r) and Vp{u — (j)) is very simple. Note that the 
random field h incorporates random torques and random 
compressions, coming from the _L and z components of 
/i, respectively. The function Vp{u — 4i{r)) represents the 
tendency of the aerogel to pin the smectic layers in a 
set of random positions <j){r), modulo the smectic layer 
spacing a, which is why Vp is periodic in its argument. 

To compute self-averaging quantities, e.g., the dis- 
order averaged free energy, it is convenient to employ 
the replica "trick" that relies on the identity logZ = 



4 



lim — — i . After replicating and integrating over the dis- 
n— >o " 

order h{r) utilizing Eq. ^2 we obtain 



^ d'rY, [k {Vlu,,)^ + B {d^u^f 
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Assuming Dq{T) is very small right at the phase 
transition the noninteracting propagator Gap {<t) = 
{ua{q)u(:j{—q)) can be easily obtained 



Gapiq) = TO {q) 5ap + AtqiG {q}' 



(13) 



withG(g) = l/ {Kq\+Bqt). 

We employ the standard momentum shell renormal- 
ization group (RG) transformation. The only novelty is 
that we will employ an infinite hyper-cylindrical Brillouin 
zone: \qi_\ < A, —oo < qz < oo, where A ~ ^ is an ultra- 
violet cutoff. We separate the displacement field into high 
and low wave vector components Uq (r) — ('^)+'"q (r)i 
where (r) has support in the hyper-cylindrical shell 
Ae^^ < q± < A, — oo < qz < oo. We then integrate out 
the high wave vector part {f), and rescale the length 
and long wavelength part of the fields with r± = r'^e^, 
z = z'e'^^, and m„ (r) = e^^u'^ (f ') so as to restore the 
UV cutoff back to A. 

Evaluating the corrections SB, SK, SA, Sg and Su and 
performing the rescalings described above, we obtain the 
following RG flow equations to one loop order: 



^=[d-5 + u; + 2x+^9,^K (15) 
^=U-3 + 3x+3253-32ff4jff, (17) 



dD{£) / 9 5, 

-dT = [d~i + u^ + ^X+YQ93-jQ9^]D 



+ -^K (.94 - gs) 



(18) 



dwje) 

d£ 



= d-5 



32 54 



3 15 
+ ( o53 - ^54 



(19) 



where 32 = A (S/A^^) 2 ^ ^ 

(5) 52, e = 5— d. These iiG flow equations have two fixed 
points: one preserving rotation invariance (53 = (74), 
which is unstable; and one with 53 = 0, gl = ^e, which 
is stable and controls the second-order phase transition. 
Analyzing the RG flows around the stable fixed point in 
the standard way leads to the critical properties, expo- 
nents, and scaling relations described earlier. 

In summary, a theory of smectic A — C phase transition 
in a field in disordered media is developed. We found 
the critical exponents to first order in the e = 5 — d 
expansion. In addition, we have made experimentally 
testable predictions for the elasticity and fluctuations of 
this system in both phases, and at the transition. 
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